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Abstract 

Evidence of how temperature takes part in the magnetic irreversibility in the intermediate state 
of a cylinder and various disks of pure type-I superconducting lead is presented. Isothermal mea- 
surements of first magnetization curves and magnetic hysteresis cycles are analyzed in a reduced 
representation that defines an equilibrium state for flux penetration in all the samples and reveals 
that flux expulsion depends on temperature in the disks but not in the cylinder. The magnetic field 
at which irreversibility sets in along the descending branch of the hysteresis cycle and the remnant 
magnetization at zero field are found to decrease with temperature in the disks. The contributions 
to irreversibility of the geometrical barrier and the energy minima associated to stress defects that 
act as pinning centers on normal-superconductor interfaces are discussed. The differences observed 
among the disks are ascribed to the diverse nature of the stress defects in each sample. The pinning 
barriers are suggested to decrease with the magnetic field to account for these results. 
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I. INTRODUCTION 



The dynamics of normal and superconductor regions in the intermediate state (IS) of 
type-I superconductors has been studied extensively by means of theoretical and experi- 
mental procedures l iol. Magnetic irreversibility has been attributed historically to the 
presence of either physical or chemical defects in the bulk of the sample and to sur- 
face and shape effects 41461] that may control both the penetration and the expulsion of the 
magnetic flux. However, this field has been recently revised on the basis of accurate dc 
magnetization measurements and high-resolution magneto-optical imaging of the flux pat- 
terns in the IS of extremely pure, defect-free type-I lead (Pb) superconductors |l3l4l5l| . The 
magnetic irreversibility has been found intrinsic for disks when the applied magnetic field 
is parallel to the revolution axis, but a fully reversible magnetization has been measured in 
these samples when the magnetic field is perpendicular to this axis or in ellipsoidal samples 
in any magnetic field configuration. The irreversibility has been thus ascribed to the pres- 



ence of a geometrical barrier t 
of the applied magnetic field 



lat depends on the shape of the sample and on the orientation 



17|. Moreover, it has been observed to be correlated with 



the formation of different flux patterns depending on which branch is explored along the 
magnetic hysteresis cycle, being bubbles for flux penetration or labyrinthine structures for 



flux expulsion |15|. However, the effect of the geometrical barrier vanishes when the 
applied magnetic field is removed and the system becomes unable to trap any flux line. 

This picture is significantly modified when the sample contains other irreversibility sources 
such as stress defects. In this situation, and only when the magnetic field is parallel to the 
revolution axis, the defects enhance the capability of the system to trap magnetic flux 
and, as a result, a large remnant magnetization can be retained at zero field Q. The 
similarity of normal-superconductor interfaces (NSI) in type-I superconductors with domain 
walls (DW) in ferromagnets led recently to suggest that the effect of the stress defects is to 
pin the NSI driving the system into a metastable state [lol. This was tested in a disk of Pb 
by studying the time-evolution of the remnant magnetization at zero magnetic field. The 
magnetic relaxation followed a logarithmic dependence which was attributed to the presence 
of a broad distribution of pinning energy barriers, in good agreement with the nature of the 
stress defects in this sample. The data obtained at low temperatures were analyzed using a 
theoretical model that describes the quantum tunneling of interfaces (QTI) mediated by the 
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formation and decay of bumps in the NSI pinned at the defects. Several physical properties 
of the system, like the average size of the energy barriers, the average characteristic lengths 
of the bumps, and the crossover temperature between thermal and non-thermal regimes of 



magnetic relaxation, were thus determined 



19|. 



However, the magnetic field dependence of the pinning energy barriers associated with the 
defects is still far from having been explored. Moreover, to our knowledge, the influence of 
temperature on either the topological or the pinning contributions to the irreversibility of the 
magnetic hysteresis cycle has not been reported yet. To investigate these effects, isothermal 
magnetization curves were measured on different samples at several temperatures and were 
analyzed in terms of reduced magnitudes defined as rrir = M/Hc and h = H/Hc, where M 
is the measured magnetization, H is the applied magnetic field, Hc(T) = Hco[l — (T/Tc)^] is 



the thermodynamic critical field 



20|, T is the temperature, and Tc is the superconducting 



transition temperature. In this representation, the resulting mr{h,T) curves of a sample 
that presents a reversible magnetization in the whole magnetic field range, like an ellipsoidal 
sample, should scale onto a single curve, rririh), whose theoretical description would only 
depend on the geometrical factor, 21|. The possible temperature dependence of the 
irreversibility of the system should be therefore reflected on the lack of scaling of the rrirlh, T) 
curves. A comparative investigation of a cylinder that only exhibits geometrical effects and 
various disks that also contain structural defects of different nature is made in this work to 
elucidate the role of temperature in the magnetic irreversibility of these samples. 



II. EXPERIMENTAL SETUP 



We have investigated four samples, labeled A to D, all made from a commercial rod of 
extremely pure lead (99.999 at.%) [22]. Samples A and B are two disks of base area 40 mm^ 
and thickness 0.2 mm that were made by cold rolling a slice of rod to the desired thickness 
and cutting out the edges to get an octagonal cross section. Sample A was further annealed 
for one hour in glycerol at 290°C in nitrogen atmosphere to reduce the mechanical stress 
associated to defects. Sample C is also a disk with a shape similar to that of samples A and 
B but it was heated above the melting point of lead (circa 327. 5°C) and cooled down to 
room temperature in a few seconds trying to induce stress defects of a nature different than 
those originated by the cold rolling protocol. Sample D is a cylinder of radius 1.5 mm and 
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length 3 mm that was produced by eroding the rod with a whetstone very slowly to avoid 
the generation of physical defects in the bulk of the sample. 

All magnetic measurements were performed in a commercial superconducting quantum 
interference device magnetometer with a low-temperature stability better than 0.01 K 23 |. 
In all cases, H was applied parallel to the axis of revolution {z direction), and isothermal 
M{H) curves were recorded at constant T values between 2.00 and 6.00 K. The samples 
were first cooled from 8.00 K (normal state) in zero magnetic field (zfc process) down to the 
desired T value, and H was then progressively increased up to an intensity slightly above 
the corresponding Hc{T) value at which M vanished, to obtain the first magnetization 
curve, Mist{H,T). Then H was swept down to —Hc(T) and up to Hc{T) again to obtain 
respectively the descending and ascending branches, M^qs{.H,T) and M.^sc{.H,T), of the 
magnetic hysteresis cycle. The number of experimental points in all the cycles was equal 
to make sure that the total time needed to record each one was always the same. Fitting 
the values of Hc{T) obtained from the Mist{H, T) curves to the theoretical expression given 
above produced values of 802 ± 2 Oe for Hco and 7.23 ± 0.02 K for Tc, in good agreement 



with standard figures that can be found in the literature for similar samples [24 1. 



III. EXPERIMENTAL RESULTS 

Panels (a), (b), and (c) in Fig. 1 present respectively the magnetization curves for samples 
A, B, and C (disks), at temperatures between 2.00 K (outermost curve) and 6.00 K (inner- 
most curve) in steps of 0.50 K, using a representation in which mr{h, T) has been normalized 
to the maximum value attained along the cycle at 2.00 K, m{h^T) = m.r(/i, T)/mr,max(2.00 
K). The three panels exhibit several common features. First, all the mist{h,T) curves 
show an initial linear regime (Meissner state) up to a certain value, h'^ ~ 0.1, at which 
the magnetization enters the intermediate state, goes through a minimum and tends pro- 
gressively to vanish at h = 1. These curves practically overlap for each sample, namely 
mist(/i, T) ~ mist{h), suggesting that the process of flux penetration is temperature inde- 
pendent to a large extent. Second, as h decreases from the normal state, a reversible region 
appears, in which all the mdes{h,T) 
called reduced irreversibility field, h* 
maximum is reached at a certain h value, below which they tend progressively to merge into 



Dranches stick together with mist{h) down to the so- 
~v\ . These branches separate at h* and stay so until a 



4 




FIG. 1: (Color online) First magnetization curves and magnetic hysteresis cycles measured on 
sample A [panel (a)], sample B [panel (b)], and sample C [panel (c)], when the magnetic field 
is applied parallel to the z direction, after a zfc process down to temperatures that increase pro- 
gressively between 2.00 K (outermost curve) and 6.00 K (innermost curve) in steps of 0.50 K (see 
legend for details). The data are plotted using the reduced m{h) representation defined in the text. 

a single curve that stretches down to h = —1. Third, a substantial amount of trapped flux 
at zero magnetic field, the so-called reduced remnant magnetization, mrem = ^dcs{h = 0), 
is observed. Fourth, as h increases back from —1, each ■m.^c{h,T) branch turns out to be 
for the most part antisymmetrical with respect to the corresponding m^csihiT) branch, and 
therefore it reproduces the behavior just described when h evolved in the opposite sense. 
In particular, each mg^c{h,T) branch merges again into mist{h) for a h value slightly larger 
than h'^, that is mi^sc{h,T) ~ mist{h) from this point on. Finally, both h* and m^era decrease 
with temperature for each sample and grow from sample A through sample C. 
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Fig. 2 shows the m{h,T) curves obtained for sample D (cyhnder) at temperatures be- 
tween 2.00 and 6.00 K in steps of 1.00 K. There are various details that distinguish this figure 
from Fig. 1. First, the shape of the hysteresis cycle changes, and in particular the Meissner 
state extends to larger h values {h'^ ~ 0.6), reflecting the fact that the demagnetizing factor 
decreases when the thickness of the sample increases (A^ ~ 9/10 for the disks 25|, ~ 1/3 
for the cylinder [2^). Second, not only the mist{h,T) curves, but the full hysteresis cycles 
obtained at different temperatures largely superimpose in the whole magnetic field range, 
to wit m{h,T) ~ m{h). Finally, h* becomes closer to 1 and mrem turns almost zero. 
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FIG. 2: (Color online) First magnetization curves and magnetic hysteresis cycles measured on 
sample D, when the magnetic field is applied parallel to the z direction, after a zfc process down 
to temperatures between 2.00 K and 6.00 K in steps of 1.00 K (see legend for details). The data 
are plotted using the reduced m{h) representation defined in the text. 



IV. DISCUSSION 

Let us start discussing the results obtained for sample D. First, the fact that rrirem tends 
to vanish indicates that the irreversible behavior in this case is only given by the presence of 
a geometrical barrier [l^ . This is reinforced by magnetic relaxation experiments performed 
at several magnetic fields along m{h) that do not show any evolution of the magnetization 
for times up to two hours, as it is expected for a defect-free sample [l^]. Each magnetic 
cycle, M{H,T), should then represent the topological hysteresis obtained in this system at 
the corresponding T value. Moreover, considering that all the cycles scale onto a single m{h) 
curve in Fig. 2, equivalent spatial distributions of flux patterns (bubbles on penetration. 
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lamellae on expulsion) should form at all temperatures for a given h value, indicating in fact 
that the topological hysteresis is thermal independent. 

Therefore, aside from the variation in the value of from one sample to another, the dif- 
ferences observed among A, B and C, and between them and D should be mostly attributed 
to the presence of structural defects of diverse nature in the disks. In these samples, as the 
magnetic field is swept, the NSI would move along, they would lose energy due to dissipa- 
tive processes, and they would be eventually pinned by the defects, leaving the system in a 
metastable state. Nonetheless, the fact that mist{h) is mostly temperature independent for 
each sample suggests that the penetration of the magnetic flux from the edges is dominated 
basically by thermodynamical processes. This means that, even if defects play an important 



role in the formation of flux patterns 



15 



271 b y shaping superconducting and normal do- 



mains via the pinning and depinning of NSI 



19| , their effect at the macroscopic scale would 



be almost negligible as the magnetic field increases from zero, indicating that mist{h) would 
represent indeed an equilibrium state for flux penetration into each disk. This is additionally 
supported by the occurrence of negligible relaxation rates when the time evolution of the 
magnetization is recorded over two hours after starting at any point along mist(/i). 

Pinning effects should become thus important in disks during the expulsion of the mag- 
netic flux along mdes{h,T) and should be therefore invoked to understand both the change 
in the area of the hysteresis cycles among samples and the thermal dependence of the cycles 
in each one. The onset of irreversibility at h* as h decreases along mdcs{h) brings about 
thermal activation processes that are needed to release the NSI from the pinning by the 
defects. In analogy with the diffusion of DW in ferromagnets [28], the probability that these 
processes occur is ruled by an Arrhenius equation, z/ = h'Qexp[—U{h)/kBT], where z/q is the 
so-called attempt frequency, which is usually of order 10'' — 10^^ Hz, and U{h) is the effective 
pinning energy barrier, which has been considered to depend actually on h. These processes 
will be relevant at the time scale of our experiments whenever the relation U{h) ~ 25A;bT 
is accomplished at a certain T value. Therefore, if this relation is to be satisfied at h*, the 
decreasing h*(T) dependence observed in Fig. 1 for all samples would imply that U{h*) and, 
in general, U (h) should be a decreasing function as well. 

The average size of the pinning barriers at zero magnetic field, Uq = f/(0), was recently 
estimated to be roughly 100 K from magnetic relaxation experiments in a system like sample 
A 19| . Similar measurements (not reported here) performed during the current investigation 
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FIG. 3: (Color online) Temperature dependence of the reduced irreversibility field [panel (a)] and 
the reduced remnant magnetization [panel (b)], normalized both to their values at 2.00 K, for 
samples A (squares), B (circles), and C (triangles). 

produced values of Uq of approximately 150 and 200 K for samples B and C, respectively. 
Thus, modeling the barrier as U{h) = Uof{h), where f{h) contains the functional decreasing 
dependence of U{h), the range of h at which the relation U{h) ~ 25A;bT holds at h* would 
be simply determined by the actual value of Uq. To be precise, as Uq increases progressively 
from A to C, h* has to increase correspondingly, as it is observed from panel (a) through 
panel (c) in Fig. 1. Furthermore, this trend is also reflected in the amount of trapped flux 
at zero magnetic fleld: as irreversibility sets in at h*, the NSI are immediately pinned by 
the defects and stay so as h decreases along mdes{h) down to zero. Therefore, the value of 
mrem wiU bc largely determined by the value of h*, and so it will also increase from A to 
C, as it is shown in Fig. 1. The different size of the pinning energy barriers can be thus 
held responsible for the change in the area of the whole hysteresis cycle from one sample to 
another. 

We can gain an insight into the thermal dependence of the cycles by looking at how h* and 
mrem chauge with temperature. Panels (a) and (b) in Fig. 3 plot respectively h*{T)/h* {2.00 
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K) and mi.eni(^)/'"^rem(2.00 K) for samples A (squares), B (circles), and C (triangles). All 
the curves in both panels show a decreasing tendency, with a slope that diminishes from 
sample A through sample C. This is a consequence of the fact that, for a fixed range of 
T values, the range of h* values that verify the relation U{h) ~ 25A;bT has to be reduced 
as Uq increases from A to C, hence the slopes of the h*{T)/h* {2.00 K) curves decrease in 
accordance. In addition, as mentioned above, the intensity of h* establishes the amount of 
mrem that can be trapped at zero field, so the relative variation of h* with T should also 
determine the range of values of mrem that can be attained at our experimental temperatures. 
Therefore, h*iT)/h*{2.00 K) and mrem(7')/'"^rem(2.00 K) will be correlated for each sample, 
as it is observed in Fig. 3. 

Finally, the reversible region that appears in the cycles at h > h* could be ascribed in this 
context to the fact that U{h) might become too small in this range of h to trap magnetic 
flux by pin ning the NSI and, as a consequence, the latter would move quasi freely through 
the sample 291]. This would be actually in agreement with previous findings of a decreasing 
magnetic-field dependence and an upper bound for the velocity at which the superconducting 
domains move in a current-induced flow experiment in a Pb slab with pinning defects 301] . 



V. CONCLUSIONS 



To summarize, we have explored the role of temperature in the magnetic irreversibility of 
a cylinder and various disks of a type-I Pb superconductor when the applied magnetic field 
is parallel to the axis of revolution of these samples. The magnetic hysteresis cycles were 
influenced by temperature in the disks, but not in the cylinder. Whereas the contribution 
of the geometrical barrier was observed to be temperature-independent for all samples, the 
pinning of the normal-superconductor interfaces in stress defects was found to be responsible 
for the different thermal dependence of the magnetic hysteresis cycle in each disk. The onset 
of irreversibility along the descending branch of the cycle was ascribed to the thermal acti- 
vation of the interfaces over the pinning energy barriers. A decreasing dependence of these 
barriers on the magnetic field was suggested to explain the fact that both the irreversibility 
field and the remnant magnetization at zero field diminished as the temperature increases. 
The variations in the values of these magnitudes and in the whole shape of the hysteresis 
cycle from one sample to another were attributed to the different size of the pinning energy 
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barriers, which in turn reflects the varied nature of the stress defects created in each sample. 

We believe that the results of this paper set the basis for a major exploration of experimen- 
tal and theoretical questions related to the magnetic irreversibility of type-I superconductors 
and, in particular, to the role of pinning of normal-superconductor interfaces in the genera- 
tion of flux patterns. For instance, the hysteresis cycle has been found to be fully reversible 
in disks when the magnetic field is applied perpendicular to the z direction fl?!, suggesting 
that pinning effects are not detected when the geometrical barrier does not play a part in 
the magnetic processes. How this barrier appears to condition the manifestation of pinning 
during the process of flux expulsion remains an open question yet. On the other hand, the 
discussion of the magnetic field dependence of the low-temperature relaxation rate of a disk 
in terms of the QTI model, assuming that the pinning energy barriers decrease with the 
magnetic field as suggested here, will be the subject of a forthcoming publication jsf] . 
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